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Abstract 



Some renormalization group approaches have been proposed during the 
last few years which are close in spirit to the Nightingale phenomenological 
procedure. In essence, by exploiting the finite size scaling hypothesis, the ap- 



o 

O . proximate critical behavior of the model on infinite lattice is obtained through 



the exact computation of some thermal quantities of the model on finite clus- 
ters. In this work some of these methods are reviewed, namely the mean field 
renormalization group, the effective field renormalization group and the finite 
size scaling renormalization group procedures. Although special emphasis is 
given to the mean field renormalization group (since it has been, up to now, 
much more applied an extended to study a wide variety of different systems) 
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a discussion of their potentialities and interrelations to other methods is also 

addressed. 
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I. INTRODUCTION 

The renorinalization group formalism introduced by Wilson in the early 70's [0,0 is by 
now one of the basic strategies to solve fundamental problems in statistical mechanics. It 
is also a very useful tool to tackle problems in several fields of theoretical physics such 
as the study of nonlinear dynamics and transitions to chaos |^, disorder surface growth 
^, earthquakes |^, among others. The conceptual foundation of the method, first laid by 
Kadanoff to qualitatively predict scaling behavior at a second-order phase transition, is 
to reduce, in a step-by-step way, the degrees of freedom of the system leaving unchanged 
the underlying physics of the problem. This reduction, carried out repeatedly through a 
renormalization recursion relation, leads the original system with a large correlation length 
(the range at which physical microscopic operators are correlated) to one with correlation 
length of unity order, where well-known methods as perturbation theory can, at least in 
principle, be used to treat the problem. Depending on the mathematical technique such 
thinning degrees of freedom can be implemented directly in the reciprocal (momentum) 
space or in the real (position) space. The former approach makes use of mathematical 
tools from quantum field theory with the crystalline system being replaced by its continuous 
limit. As a result, the so-called e— expansion proposed by Wilson and Fisher ^ (and further 
developed by using techniques of renormalized perturbation theory (§0) provides analytical 
and quite well controlled asymptotically exact results for critical exponents (despite being 
unable in predicting values of critical points, critical lines and phase diagrams). On the other 
hand, the more intuitive real space version of the renormalization group works directly in the 



position space. It was introduced by Niemeijer and van Leeuwen [|lOl and several different 



techniques have been proposed and applied to a great variety of statistical models |]ll 
The real space renormalization group (RSRG) has since become an important apparatus in 
studying critical phenomena. 

The main feature of the renormalization group is to obtain, from the renormalization 
recursion relations, flow diagrams in the parameter space of the system, i.e., the space 



spanned by the different kinds of interactions present in the Hamiltonian. The attractors 
(trivial fixed points) and their basins in the flow diagram correspond to regions of different 
thermodynamic phases. These regions are separated by critical frontiers associated to semi- 
stable attractors (relevant fixed points) which determine the universality class of its critical 
exponents. More unstable attractors located on the critical surfaces lead to multicritical 
frontiers (or, as it happens to be more usual, just multicritical points). There can also be 
a first-order frontier line linking a multicritical point to a first-order fixed point located at 
infinite values of the Hamiltonian parameters. Moreover, the singularities of the critical 
frontiers on the multicritical regions are characterized by the crossover exponents. So, from 
the above general view, the topology of flux diagrams allows one to achieve qualitative 
information about the critical behavior of the system, e.g., universality, order of transitions, 
crossover and multicriticality, among others. Regarding the real space renormalization group 
approach, which will be the subject of this work, one can also say that in most cases rather 
accurate quantitative values for the exponents and transition lines are obtained. 

Just for completeness we should say that, from a different and more modern point of 
view, the renormalization group is essentially based on the fractal structure exhibited by the 
system at the critical point. A fractal, named by Mandelbrot [|r^], describes structures con- 
sisting of parts of any size (i.e., in any scale) which are similar to the whole. Consequently, 
the renormalization group commonly establishes a recursion relation transforming the pa- 
rameter space of the Hamiltonian into itself, leaving thus the underlying physics unchanged. 
Due to this continuous scale invariance a fractional dimension is ascribed to the fractal (and, 
as a consequence, to the critical exponents, as we shall shortly see in this Introduction), in 
contrast to the integer dimension of the translationally invariant Euclidian space. There 
are, however, systems exhibiting discrete scale invariance (a weaker form of scale invariance 
symmetry which is not continuous) where dimensions or exponents are complex. Even in 
such cases the renormalization group formalism turns out to be quite useful 0. The in- 
terested reader is referred to the recent review by Sornette for further details concerning 
discrete scale invariance and complex dimensions |jl3| . 



The renormalization group procedure has been extensively reviewed in hterature (see, 
for instance, the excellent reviews in References |0,0,^|lT|; see also a more recent text by 
Yeomans on this subject |1^). Regarding the RSRG several techniques have been proposed 



and applied with success to various problems. Among the simplest ones we have decimation 
I3| , Migdal-Kadanoff |]TB[ and Niemeijer-van Leeuwen cells |]TD[ approaches, which have 



also been discussed in the reviews referenced above. More accurate procedures include 
Monte Carlo renormalization group [|l7l, correlation function preserving renormalization 



group ||T8| and phenomenological renormalization group ||19|- It is not the scope of the 
present work to discuss all the methods above, rather just the phenomenological ones will 
be treated. Nevertheless not less excellent reviews on Monte Carlo renormalization group 



| 20| and correlation function phenomenological renormalization group |^ are addressed to 



the reader. 

Despite its success, the ordinary RSRG has some drawbacks in the core of its imple- 
mentation. Although exact results are in general achieved in classical models for dimension 
d = 1, in most cases ioi d> 1 some uncontrollable, and rather obscure, approximations have 
to be done in order to properly obtain the renormalization recursion relations. In addition, 
one has also to deal with the awkward "proliferation of parameters" , usually encountered in 



approaches such as decimation |15 



One class of RSRG methods, free from the drawbacks presented above, is the so-called 
"phenomenological renormalization" or finite size scaling approach. In these methods one 
computes (exactly) thermodynamic quantities Pl and Pli for two different finite systems 
and, from the expected scaling relation obeyed for this quantity P in the limit where L'{< L) 
tends to infinity, renormalization group recursion equations can be obtained. In this way the 
critical properties of the infinite system are obtained (approximately) from the knowledge 
of the corresponding properties of its finite lattice counterpart. 

The finite systems used in these phenomenological renormalizations depend strongly on 
the geometry of the lattice. For hypercubic lattices there are two specific geometries of 
particular interest: 
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i) a finite system in all directions consisting of a hypercube of side L in a d-dimensional 
lattice; 

ii) a system infinite in one direction and finite in the other d — 1 dimensions with cross- 
sectional area characterized by the length L. 

Other geometries, such as systems infinite in two dimensions, can also be considered. 
However, due to the difficulty in obtaining exact solutions even for two-dimensional models, 
such geometry will not be discussed here. On the other hand, boundary conditions applied 
to the finite systems depend on the particular renormalization group strategy. The common 
ones are the periodic boundary conditions and, in this case, just bulk criticality is studied. 
There are, however, some methods which allow to study bulk and surface criticality by taking 
free boundary conditions in the finite lattices. Sometimes helical boundary conditions are 
imposed in geometries ii) in order to obtain a sparse transfer submatrix leaving the problem 
numerically more tractable p2| . 



To illustrate, in a quite simple way, the procedure of writing down a RSRG recursion 
relation for finite size systems let us take, without loss of generality, the magnetic language 
by means of the spin-S Ising Hamiltonian model defined as 

JV 

n = -jY. s,s,-hY,s,, (1) 

<i,j> 1=1 

where the first term is a sum over nearest neighbors pairs of spins < i,j > (which might 
be conveniently extended to more distant neighbors) located on an arbitrary lattice and the 
second term is a single sum over all the A/" ^ oo spins of the lattice. J is the exchange 
interaction which is ferromagnetic for J > or antiferromagnetic for J < and h is an 
applied external magnetic field. The spin variables Si take values Si = S,S — 1,...,—S 
where S* is a fixed integer or semi-integer positive number. This model exhibits a second- 
order phase transition ai K = Kc and H = 0, where K = (3J is the reduced exchange 
interaction, H = (3h is the reduced external magnetic field, (3 = l/ksT and fc^ is the 
Boltzmann constant. At zero external field the correlation length ^ close to the transition 
is expected to diverge as 



e^ier% (2) 

where e = K — K^ ^ and z/ is the corresponding critical exponent. 

From the Kadanoff block construction it is shown that close to the transition the sin- 
gular part of the free energy per site /s(e, H) is a homogeneous function of its thermodynamic 
variables (as stated by Widom in 1965 |2^) i.e., 

fsie,H)=r'fs{f/''e,iy-H), (3) 

where i is an arbitrary scaling factor, d is the dimensionality of the corresponding lattice, 
and yn is the magnetic critical exponent related to the magnetization m through m = H^'^, 
with 6 = yu/i)- — Vh)- Now, any other quantity P, obtained as derivatives of the above free 
energy, will behave as power laws 

P ^ |e|-" , (4) 

where a is the critical exponent of the quantity P. For example, the magnetization 

Tr exp(—p7t) 
is obtained, close to criticality, from 

m{e, H) = -^^%-^ = r'^+^«m(£i/'^e, iy^H) . (6) 

oH 

At H = and choosing i~^/'^ = e (e > 0) Eq. (H) gives, after comparing to Eq. (||), 
a = P = —{d — yH)t^- Similarly, for the zero field specific heat C = — T {d"^ fs/dT'^)^^^ and 
zero field magnetic susceptibility x = {9^^/dH)jj^Q it is easy to obtain a = a = 2 — du, and 
(J = ■y = {2yH — d)^, respectively. These power laws are a signature of the continuous scale 
invariance of the system at criticality. 

On the other hand, a finite system of linear dimension L close to K^ will of course present 
a non-diverging characteristic length and its correlation length will, at the maximum, be 
given by ,^ ~ L. Inserting this relation in Eq. (0) and substituting for e in Eq. (^) one gets 
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p ^ L"/^ = L^ , (7) 

where = cr/z/ is called the anomalous dimension of the quantity P. As is commonly 
a fractional number it is clear again from the above relation the effects of continuous scale 
invariance on dimensions and exponents of critical systems. 

According to the finite size scaling hypothesis, the generalized scaling relation obeyed by 
any thermodynamic quantity P taking into account the finite size L of the system can be 
expressed as p4| , p5| 

P(e, H, L) = e'^Pif/^e, &"H, i'^L) , (8) 

where L — > cxo is the linear dimension of the finite system and e ~ 0, if ~ 0. By desiring 
the right hand side of Eq. (^ to be the respective quantity P for a smaller system of linear 
dimension L' we choose for the scaling factor £~^L = L' or 

Expressing now the variables i^^'^e = e' = K' — K^ and &" H = H', Eq. (W) can be rewritten 



as 



PLiK',H') _ PLiK,H) 

where the system size is now given as a subscript instead of a variable (to be more apparent) 
and primed quantities refer to the smaller system V . This is the basic equation from which 
a mapping {K, H) — > {K' , H') is obtained with the rescahng factor defined in Eq. (H). It is 
clear that the only source of inaccuracy in this relation resides in the finite values taken for 
L and L' in order to get exact calculations for Pl and Pl>. Nevertheless, it has been noted 
that such approaches are applicable to rather small systems of no more than a few degrees 
of freedom. Of course, the bigger the finite systems, the better the results achieved, and it 
turns out that sometimes they exceed in accuracy more conventional methods such as series 
expansion, e-expansion or other renormalization group procedures. 
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Different quantities P will of course generate different approaches. The first procedure, 
and in our opinion one with the most accurate results, was proposed by Nightingale |T^ by 
using the correlation length computed for finite lattices consisting of infinite strips with finite 
width, i.e., the geometry ii) discussed above. From the basic ideas of the renormalization 
group theory it is well known that the transformed system with less degrees of freedom has 
a smaller correlation length which is the old one reduced by the rescaling factor i, namely 

e{K',H') = ^-^^. (11) 



From equations (^ and ([TT|) we then get 

ejK'.H') ^ i{K,H) ^ ^^2) 

Ij Ij 

This means that for the correlation length the anomalous exponent is equal to unity and the 
above equation was interpreted by Nightingale as a renormalization group transformation of 
the infinite system from which fixed points, critical exponents, etc., are obtained. Since Eq. 
(P^) involves just a single function, we obtain just one scalar recursion relation instead of the 
usual multidimensional recursion relation on a set of coupling constants such as {K, H) in 
the present example. So, flow diagrams analyses are not possible when the Hamiltonian has 
more than two parameters. However, second-order transition lines and estimates of critical 
exponents can be obtained by determining a mapping 

K^K' = R{K;r) (13) 

for a fixed value of r = H/K, where the recursion relation R{K; r) is, in principle, obtained 



from Eq. (|T3). The critical line is given by the fixed point solutions K' = K = K* = 
R{K*]r) as a function of r and estimates of the critical exponent is obtained from the 
linearization of R around K* 

K' -K* = Xt{K - K*) , (14) 

where the thermal eigenvalue At is given by 
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dK' 



dR{K; r) 



dK 



dK 
with the corresponding thermal critical exponent 



(15) 

K* 



f/" = Xt . (16) 



Estimates of the magnetic critical exponent are extracted from 

dH' 



PJH 



dH 
from which, using Eg. (p!2D , gives 



(17) 

K* 



The second derivatives appear in obtaining yn since ^ is commonly an even function of H. 
The procedure outlined above is quite general when using only one recursion relation and 
should be applied to any function R{K\ r) (coming, for example, from other renormalization 
group scheme) as well as to any couplings (not necessarily an external field). 

Other thermodynamic functions (such as, for instance, specific heat or magnetic suscep- 
tibility) can also be used to study the critical properties of statistical mechanical systems 
through Eq. (pUj). The main problem in these cases is that the exponent is in general not 
known. The method can, nevertheless, be implemented by force through the use of three 
different finite systems (L, L', L") and taking the estimate of (j) the value that yields the 
same fixed point solution for the recursion relation (|1^) from (L, L') and (L', L") clusters, 
respectively p6[. In the present work we will be concerned only with methods where has 



a known value. 

By taking P as the order parameter of the system and together with mean field calcula- 
tions one can obtain the so-called mean field renormalization group (MFRG) approach ||27|| . 
Here, it is possible to overcome the difficulty of not knowing the order parameter anomalous 
dimension and finding out a recursion relation free from any exponent. This procedure 
will be discussed in more detail in the next section. It is also possible, in some models hav- 
ing more than one order parameter, to obtain complete fiux diagrams in the Hamiltonian 

11 



coupling space due to additional recursion relations. Moreover, by taking three clusters at 
the same time, it is allowed to study bulk and surface critical behavior since the employed 
finite lattices must have open boundary conditions. 

It is clear from Eq. ( p!OD that in carrying out phenomenological renormalization group 
calculations one desires to choose quantities P in which their anomalous dimension are 
known a priori. A more recent proposal, for particular quantities having = 0, has been 
introduced by de Oliveira in the finite size scaling renormalization group (FSSRG) method 
p8|^0|. This method is also capable to circumvent, in some cases, the lacking of fiow 
diagram lines in the parameter space by considering more than one quantity at the same 
time and obtaining additional recursion relations, as in some versions of the previous MFRG 
approach, even for models presenting just one order parameter. This approach will also be 
discussed later on. 

The purpose of this work is thus to review some phenomenological real space renormal- 
ization group procedures namely, the mean field renormalization group, the effective field 
renormalization group (EFRG) and the finite size scaling renormalization group. We will 
not address here further details on the Nightingale phenomenological renormalization since 
it has already been widely discussed in literature (see, for instance. Reference [^ and refer- 
ences therein). Special emphasis will then be given to the MFRG since it has been extended 
and improved in different ways and applied to a vast class of models (classical, quantum, 
geometrical, pure and diluted, static and dynamic). We will also not discuss here exten- 
sively the results achieved in each of its applications. We will be mainly concerned here in 
all the versions of the method and its general formulations, some of them not previously 
published in the literature. Just a simple list of references in alphabetical order (without 
any comments) on works involving the use of MFRG up to the year 1993 has already been 



given by Croes and Indekeu [32 



The review is arranged as follows. In the next section we present the formalism of the 
mean field renormalization group for the study in static and dynamic critical phenomena. 
In Section III the exact results from the MFRG are presented. Section IV gives an up to 
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date application of the method in static problems, and Section V does the same concerning 
the dynamic critical behavior. Related phenomenological renormalization group approaches 
are presented in section VI and some final remarks are discussed in section VII. 

II. PHENOMENOLOGICAL MEAN FIELD RENORMALIZATION GROUP 

APPROACH 

The phenomenological mean field renormalization group approach uses the order param- 
eter as the quantity to be renormalized. In magnetic models it is given by the magnetization 
of the system which can generally be written as m(K, H), where the vector K represents all 
the reduced coupling constant interactions of the Hamiltonian and H is the reduced external 
magnetic field. It means that the first sum in Eq. (|l]) may now not be restricted to first 
neighbors but may also include more distant interactions as well as additional fields. For in- 
stance, when K is written as K = {Ki, K2, K3, K4^, ...), Ki should represent nearest neighbor 
interaction, K2 next nearest neighbor interaction, K3 a crystal field interaction, K4 a four 
spin interaction, and so on. Since the MFRG has been proposed by Indekeu, Maritan and 
Stella in 1982 [^ the method has been widely used in treating several statistical mechanical 



problems as well as frequently improved in different aspects of its formulation. This section 
is devoted to review its basic assumptions and discuss some new implementations. 

A. MFRG 

In the original MFRG one first considers two clusters of interacting spins containing A^ 
and N' sites, with A^' < A^. The surrounding spins of these clusters are fixed to a value b and 
b', respectively, which, in a sense, can be viewed as effective magnetizations representing the 
remaining spins of the infinite lattice and behave as a symmetry breaking field in each cluster. 
The magnetizations per site mAr(K, H; b) and mAr/(K', H'; b') can exactly be computed from 

"^-(^' ^' '^ = TreM-mN) ' ^''^ 
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where the trace is taken over the ensemble defined by the cluster Hamiltonian 7ijv(K, H] b) 
and Sj is the corresponding spin operator. Similar expression holds for the smaller cluster 
A^'. The usual mean field type approximation is obtained by assuming m7v(K, H;b) = b (or 
mN'(K.',H';b') = b') from which m = b (or m' = b') is solved self-consistently at i7 = 
{H' = 0) and the criticality is identified with the bifurcation point in the above equations. 
Although this mean field approximation improves the critical couplings as the cluster size 
gets larger, the critical exponents are always the classical ones (in particular, for the Ising 
model described in Eq. (|l|) one has a = 0, [3 = 1/2 and 7 = 1 for any lattice dimension). 

Instead of the naive assumption of just equalling the symmetry breaking fields to the 
magnetization per spin for each cluster, the MFRG assumes that these magnetizations are 
related through the scaling relation given by Eq. (|^) for the order parameter of the infinite 
system close to criticality ||10[ 

(20) 



mjv'(K', H'; b') = £'^-^^miv(K, H; b) , 



where the rescaling factor (. can be defined by the usual way 

1={N/N'f'\ 



(21) 



which is equivalent to definition (^ for symmetric blocks of size N = L'^. As Eq. (^) should 
hold for systems near the critical point, the magnetizations for each finite cluster must be 
very small. Such condition is achieved by letting 6 <^ 1 and 6' <^ 1, as well as H <^ 1 and 
H' <^ 1, since any finite system has no spontaneous magnetization. Equation (|20|) can thus 
be expanded and, to lowest order in b, b' and H, H' gives 



gN'{'i^)H' + /jv'(K')&' = ^"-^^/^(K)6 + £'^-^«^jv(K)i7 , 



(22) 



/^'(K') = 


dniN,{V.',H'-b') 
db' 


H'=0,b'=0 


^7V'(K') = 


dmNiK',H'-b') 
dH' 


H'=Q,b'=0 



/^(K) 



9n{K) 



dmN(K,H;b) 



db 



dmN(K,H;b) 
dH 



H=0,b=0 



H=0,b=0 



(23) 



(24) 
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As b and b' are also viewed as magnetizations and are very small they are assumed to 
satisfy the same scaling relation given by Eq. (|20|), i.e., 



b' = i^-y^b . (25) 



From equations (^) and (pSD one then gets 

/^.(K') = /^(K), (26) 

obtaining by equalling the coefficients of b' and 

gj,,{K')H' = i^~y"g^{K)H , (27) 

obtaining by equalling the external field terms. Equation (|26|) , which is independent of any 
exponent, is interpreted as a recursion relation for the coupling constants K according to the 
ordinary MFRG, while equation (^7]) , as shown below, is used to estimate the magnetic crit- 
ical exponent yn (the recursion relation does not depend on the external field). Clearly, for 
a multidimensional Hamiltonian parameter space this scalar recursion relation (p6D does not 
provide complete flux diagrams. This sole equation can, however, be used to estimate critical 
exponents and locate critical surfaces by considering a mapping Ki ^ K[ = R{Ki, ri, r2, ...) 
for fixed values of the ratio ri = K2/K1, r2 = K^/Ki, ..., where K = {Ki,K2, ...) and the 
function R is, in principle, obtained from the rescaling relation (|26|). By computing the 
fixed point solution Ki = K[ = Kl = R{Kl,ri, r2, ...) as a function of r one determines the 
critical surface, represented here as K* = {Kl,ri,r2, ...), and the corresponding exponents 
are achieved through 

Ini 



V 



In At 
as stated by Eq. (|16D for the correlation length critical exponent where 



(28) 



A -^^^ - r^^{^^^'\ 1 roq^ 



and 
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gN'{K'')=i''-^y"gr,iK*) (30) 

for the magnetic critical exponent yn, where the relation H' = i^^H has been used in Eq. 

This approach has been applied to a great variety of statistical models and quite good 
results have been obtained even by taking the simplest choice for the clusters, namely, 
A^' = 1 and N = 2. Some exact results are also obtained from this method when treating 
some specific systems. A summarized discussion of such studies will be given in the next 
sections. However, in order to taste the simplicity of the method and its potentialities let 
us apply it to three different problems by taking the smallest clusters depicted in Fig. 1. 

spin- 1/2 nearest neighbors Ising Model 

The Hamiltonians Tii and 7^2 for the clusters of one and two spins shown in Figs. 1(a) 
and 1(b) can be written as 

Hi = -cJ'h'Si - h'Si , (31) 

^2 = -JS1S2 - (c - l)Jh{Si + ^2) - KSi + ^2) , (32) 

where Si = ±1 and c is the coordination number of the lattice. From the definition of the 



magnetization (|1^) one easily has 

mi = tanh{cK'b' + H') , (33) 

smh.[2{c - l)Kh + 2H] 
""' ~ cosh[2(c- 1)7^6 + 2iJ] + e-2i^ ' ^ ^ 

which, for small values of b,b',H and H' reduce to 

mi = cK'b' + H' , (35) 

2(c-l)Kb 2H 

The usual mean field approximation assumes mi = b' and m2 = b yielding, at zero external 
field, Kc = 1/c and 1 + e~'^^'= = 2(c — l)Kc. One can note that from both clusters there is a 
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(wrong) finite critical temperature Tc ^ for tlie one-dimensional model (c = 2). For a future 
comparison to the two-dimensional model (c = 4) the later expressions give Kc = 0.250 and 
Kc = 0.286 for A^' = 1 and N = 2, respectively. Moreover, by further expanding equations 
(p3|) and (0) up to order b' and b^ it is easy to show that the corresponding critical 
exponents take the classical values, independent of the dimensionality of the lattice. 
Reminding now equations (|2^^7D the corresponding MFRG recursion relation 



gives 



--^^^ (3. 



A-. = iln^. (38) 



where Kc = K' = K is the non-trivial fixed point solution of the equation (|37|). The 



above result is the same as that obtained from the Bethe approximation [33|. The critical 



exponents according to equations (E9|) and (|30|) are 



£i/- = l + i^Jlln^, id-2y, = i^Jl £ = 2'/''. (39) 

2(c-l) c-2 ' c ' ^ ^ 

In this case, for the one-dimensional model (c = 2) one has the exact critical temperature 
Tc = and critical exponent yn = ^ (these values are in fact obtained for any cluster sizes), 
while the critical exponent z/ = 1 is only achieved by comparing two chains with A^ and A^ — 1 
spins in the limit A^ -^ oo. Numerical values for the two-dimensional model are Kc = 0.347, 
u = 1.667 and i/h = 1.415, which are, when compared to the previous mean field results 
(and without further effort), closer to the exact ones, mainly for the critical coupling (see 
Table D . The corresponding results for the d = 3 lattice using the present clusters are listed 
in Table 0. 

Site percolation 

The site percolation |^ is a geometrical problem defined on a regular infinite lattice 
where each site can independently be occupied or empty with probabilities p and 1 — p, 
respectively (p is also viewed as the concentration of the occupied sites). A cluster is formed 
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by grouping together nearest neighbors occupied sites. At p = 1 one has an infinite cluster 
finking all sites. On the other hand, for small p one has an infinite number of independent 
finite clusters. Then there is a critical concentration pc above which the size of the clusters 
(at least one) becomes infinite. The mean linear size of the clusters is the correlation length 
^ which, close to the critical threshold diverges as 

^^Ip-Pc]-"" , (40) 

where Up is the percolation correlation length critical exponent, pc depends on the topology 
of the lattice and Up depends only on the lattice dimensionality. 



Let us now consider the same clusters (a) and (b) in Fig. 1 for this problem [^3| , p| . The 
order parameter P^ is defined by the probability of a particular site to belong to the infinite 
cluster, b and b' will be now the effective values of the probability, for surrounding sites of 
each block, to belong to the infinite cluster. For the one site cluster (A^' = 1) its probability 
of belonging to the infinite cluster is p'b' times the number c of surrounding sites, i.e. 

Pi = cp'b' . (41) 

In a similar way, the probability of any of the two sites of the N = 2 cluster is, up to 
first -order in b, given by 

P^ = p{l - p){c - l)b + 2p'^{c - l)b = p{l + p){c - 1)6 . (42) 

The first term on the left hand side of the above equation takes into account the contribution 
of the configuration where one site of the cluster is occupied and the other one is empty, and 
the second term the contribution when both sites are occupied. In all terms just configura- 
tions with one surrounding site occupied with probability b are taken since configurations 
presenting two or more mean field sites will contribute to the order parameter a factor of 
order b^ or greater. Mean field approximation (Pi = b' and P2 = b) gives pc = 1/c for the 
smaller cluster and Pc = ^ { y fzf ~ 1 ) ^01 the bigger one, while the MFRG furnishes the 
recursion relation 



p'c = p{l+p){c-l) , (43) 



from which one obtains the non trivial fixed point p = p' = Pc 



Pc = ^ , (44) 



again identical to the Bethe approximation |^. The corresponding critical exponent is 
given by 

£i/-P = £±1 . (45) 

c 

As in the previous case, for the one-dimensional model the mean field approximation results 
Pc < 1, contrary to the exact value pc = 1. This exact result is however obtained from the 



MFRG recursion relation (|4^). The correct value oi Up = 1 is achieved only in the limit 
of infinite chains. For the two-dimensional square lattice equations (|4^) and (^) furnish 
Pc = 1/3 and u = 1.55, which should be compared to the series values pc = 0.59 and 
u = 1.33 1^. The results are analogous for d = 3. In addition, the related problem of 
bond percolation where bonds are independently present linking two neighboring sites with 
probability p or absent with probability 1 — p can similarly be treated from this formalism. 
It turns out that the results are identical to those above for the site problem when using the 
smallest clusters. Only by taking bigger finite systems the approach is able to distinguish 
between the site percolation and the bond percolation problems. 

Directed self-avoiding random walk 

Let us briefiy mention the directed self-avoiding random walk since it is, in some sense, 
related to the site percolation problem treated above. In order to define the directed self- 
avoiding random walk (DSARW) let us consider a c?- dimensional hypercubic lattice |[47|| . 
The DSARW is the geometrical object built by starting at a given point (the origin) and 
advancing at random through steps linking nearest neighbors sites, with the constraint that 
any step be only along the positive directions of the cartesian axis. Letting w be the fugacity 
associated with a step, one has that for w larger than a critical value Wc the walk will, on 
average, extend to an arbitrarily large distance. In this case one can distinguish two different 
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directions: one parallel (||) defined along the bisectrix of the angle formed by the cartesian 
axis, and another perpendicular (±) taken orthogonally to the previous one. It turns out that 
as w approaches Wc the root mean square of the parallel and perpendicular displacements 
diverge as power laws of the form {w — WcY'<^ and {w — WcY-^, respectively, with z/y and u^ 
being the corresponding critical indices. Figs. 1(c) and 1(d) illustrates the smallest possible 
clusters which can be used by the MFRG in order to treat the DSARW on a two-dimensional 
lattice. 

Defining now the order parameter P/v as the generating function of all walks that begin 
at the origin site and end at the mean field sites (with b being the effective value of this 
order parameter for the surrounding sites) one arrives at 

Pi = 2w'b' , (46) 

since there are two ways with probability w'b' of linking the origin of this smallest cluster 
to its two boundary sites and 

P2 = 4:w% , (47) 

where for this cluster one has four different ways of linking the origin to the three boundary 
sites all having the same probability w'^b. The MFRG recursion relation is then w' = 2w^ 
from which one has the exact non trivial fixed point w^ = 1/2. Exact values for the critical 
indexes (which are more laborious to be computed) z/|| = 1 and i'±_ = 1/2 are also obtained 



(in fact, for any cluster size, as well as for any dimension d. See Ref. [Q for more details). 
Although in all models above the calculation of the order parameters could be readily 
performed for the smallest clusters, larger systems will drastically increase the analytical 
and computational effort. Tables | and ^ contain the results obtained for the spin-1/2 Ising 
model with nearest neighbor interactions up to the greatest size treated so far (L = 6 for 
the two-dimensional model and L = 3 for the three-dimensional one), and Table fT| the 
corresponding values for the two-dimensional site percolation problem. A more detailed 
discussion of the values in these tables will be given in the next subsection. 
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Despite its simplicity the usual MFRG treated above deserves some critics regarding the 
use of larger clusters. Since the symmetry breaking fields b and b', in some sense, take into 
account the remaining of the infinite lattice, the scaling imposed for the magnetizations rriN' 
and rriN is the relation obeyed by the infinite system order parameter (note that equation 
(PPP comes from the Wilson renormalization group strategy). On the other hand, the final 
expression of the MFRG, namely Eq. (|26|) , is independent of the symmetry breaking fields. 
From this equation the functions f^ and /n' can be thought of being quantities satisfying 
the finite size scaling relation (|10D with anomalous dimension = 0. However, for very 
large clusters such relation for f^ and J'n' (with = 0) is not true. This is most easily 
seen by taking the limit of A^ and A^' going to infinity. It is clear that while the computed 
magnetizations are bulk quantities, the surrounding effective magnetizations on each cluster 
behave as surface quantities and must properly scale as surface fields. In this limit Eq. ( pSj) 
is no longer valid once yn is the bulk magnetic critical exponent. As a consequence, the 
procedure will not reproduce, in general, the exact results as the size of the systems tend to 
infinity. 

B. Surface-Bulk MFRG (SBMFRG) 



A way to handle this situation has also been proposed by Indekeu and co-workers p9[ 



For large clusters, except for spins at corners or edges, almost all boundary spins are subject 
to the same effective field which is proportional to the surrounding fixed magnetizations. 
According to a rigorous finite size scaling derivation these surrounding fields do not scale 



like a magnetization itself but like a surface field |36|. So, for b and b' one has 



b' = P'^^b , (48) 

where yus is the corresponding surface field critical exponent. One can note that now the 
factor £^~2/H ig Yio more eliminated from equations ( ^2]) and ( ^HD and a simple formulation, as 
expressed by Eq. (P^D, is not available. It is needed then three clusters A^, A^' and A^" (with 
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N > N' > N"), instead of two, in order to get a renormalization group transformation. 
In this case, from mAr(K, H; b), mN'{K', H'; h') and mjv//(K", H"; h") one computes /jv(K), 
/^/(K'), /7V"(K") and from Eq. (mj aX H = H' = H" = together with 



b' = i%"^,b , b" = i%%„b' , (49) 



it is easy to obtain 



hiK') = e%-j^r'"' MK) , (50) 

f^4K")=ip^^r'"'f^,{K'), (51) 



where 



iNN' = {N/N'f" , l^.^n = {N'/N")''" . (52) 



Following the well-established optimization strategy ||2^ the exponent d — yu— yus is deter- 
mined self-consistently by imposing that the two different mappings with rescaling factors 
Cnn' and dN'N" possess the same fixed point K = K' = K" = K*. This procedure allows 
one to obtain unambiguously the index d — yu — yns while two different estimates for the 
critical exponents v, yn and yns are achieved from each of the renormalization group trans- 
formations ( ^OD and (|5l|) (there are in fact three different estimates since from equations 
(0) and (|5TD one can construct an additional relation between /^v" and /at). 

At this point it is worthwhile to see a quantitative and comparative application of these 
methods in the well known control system square two-dimensional Ising model. Table | 
shows the critical temperature and critical exponents for different growing finite size square 
clusters. It is clear that a better accuracy is achieved from SBMFRG as compared to the 
usual MFRG, mainly for larger clusters, where the approach to the exact results is rather 
apparent. However, concerning small clusters, the distinction of both methods is not so 
profound. First, in assumption (EBI) the ordinary MFRG assumes a zero value for the 



critical index d — yu — yns (the exact one being —0.375). From the SBMFRG one gets 
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the best value d — yu — Vhs = —0.217 for the biggest systems while for the simplest choice 



A^ = 4, N' = 2 and A^" = 1 one has —0.124 E^], which is a quite small value. Second, for 



small clusters there is no even a clear distinction between surface and bulk quantities. This 
means that the fact of the SBMFRG being exact when A^, A^' and A^" -^ oo does not make 
it more correct than the original MFRG for small systems. Which method should be applied 
to a particular problem must be decided on practical grounds. Table || lists the results for 



the d = 3 Ising model and Table |T| the values for the d = 2 site percolation problem. 

Another point deserving some discussion concerns the definition of the rescaling factor 
i. As pointed out by Slotte |5^ the rescaling factors given by equations (|2l|) and (|52D 



are, to some extent, arbitrary. Different definitions strongly affect the results of the critical 
exponents for small clusters. Slotte suggested a definition of the rescaling factor where the 
length of the cluster is measured by the number of interactions, including those with the 
surrounding mean field sites, instead of just the effective number of sites, i.e., 

L-(^J:L-^'. (53) 

where the sum runs over cartesian directions and Lj are the number of bonds along each 



direction. The one spin cluster has then L = 2 and the two-spin cluster L = 6Jd/{9d — 5). 



Such definition improves substantially the critical exponents |£0[- There are, however, in the 
original MFRG, some quantities that are independent of i. One of them is clearly the fixed 
point K*, as can be seen from Eq. (pB]). The other one is the product lyi/H- By abandoning 



the definition given by Eq. (|2T|) the relation (|30[) can be rewritten as 

N'gNiK*)=r'y"Ngr,{K*) , (54) 

which gives 

yn = 2 ^^ > (55) 

where £'^ in ( ppp has been replaced by N/N'. While the exponents z/ and yn separately 
depend on the rescaling factor, it is easy to see now from equations (|28|) and (|55D that the 
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product vyn is independent of i, removing thus a principal weakness of the MFRG. In the 
unified SBMFRG this is no longer strictly true, but should still hold to a first approximation. 
It is then expected that the critical exponent product vyn will present a better precision 
than the exponents themselves. Indeed, from the last column of Table I one can see that 
the best product vyH from SBMFRG is within 1.8% of error from its exact value, whereas 
the error in critical exponents ranges between 7.8% — 9.4%. The figures for MFRG are, 
respectively, 3% and 14% — 16%. Roughly the same behavior is noticed in the data of Table 
m for the (i = 3 Ising model. 

Similar results are also obtained in an extended version of the present approach by 
including surface and corner critical exponents in two-dimensional models |51||. In such 
d = 2 systems the distinction between surface and corner fields is easily made. All one has 
to do is to include an additional corner field critical exponent ync in the formalism above. 
Table |I^ quotes the results from Croes and Indekeu [^ for the two-dimensional square 
Ising model. In this reference one also finds critical values for other two-dimensional lattices 
namely, triangular, honeycomb, hexagonal as well as the square lattice with next nearest 
neighbors interactions. On the other hand, additional edge fields will be present in three- 
dimensional finite lattices. However, bigger finite systems will require much computational 
effort and inclusion of edge fields would not be so relevant in rather small systems. 

C. Extended MFRG (EMFRG) 

In the previous subsections just Hamiltonian models with one parameter have been 
discussed where the trivial one-dimensional renormalization group flux is obtained. Let us 
now turn our attention to the case where the Hamiltonian presents more than one parameter 
(or coupling constant). It is clear that in this case the complete renormalization flow diagram 
is not defined through just the one-dimensional recursion relation as given by Eq. (p6[). 
In the particular two-dimensional case where we have K = {Ki,K2), which is the most 
studied in literature and will be treated in this section, estimates of critical lines and critical 
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exponents for models presenting only one order-parameter are often obtained through the 
mapping Ki —>■ K[ = R{Ki, r), for fixed values of the ratio r = K2/K1, as has already been 
discussed in the Introduction. That is what one finds, for instance, in the random Ising 
model or the transverse Ising model. However, depending on the system we are studying, 
other additional patterns may occur where the above simple scheme cannot be applied for 
a two- fold coupling space (in some models even for a one order-parameter system). In these 
cases the thermodynamic quantities should be considered in a proper way. This subsection 
is devoted to present a detailed discussion of what has been done in overcoming such debility 
of the method as well as to present an attempt to unify the approaches done so far by using 
the MFRG. The formalism is intended to be quite general. However, whenever possible, 
some reference to particular models (without presenting results; the corresponding reference 
will be listed in section IV) exhibiting the behavior under consideration will be made just 
for clarity and as a matter of example. 

1. one order-parameter 

For systems like, for example, the antiferromagnetic spin-1/2 Ising model in an external 
field (recall Eq. (1)) 

n = JY.s^SJ-hY,s,, Si = ±i, (56) 

where Ki = f3J and K2 = jSh and the spin-1 Blume-Capel []53| , |5^ model 



n = -j Y. SiS,-Aj2s^ , s, = ±i,o, (57) 

where Ki = j3J and K2 = (3 A (the latter interaction being the crystal field anisotropy), 
one has two couplings K = {Ki,K2) and just one order-parameter, namely, the difference 
m_|_ — m_ of the sublattice magnetizations m^ and m_ in the former case, and the mean 
value of the spin < Sj > at any site of the lattice in the later model. Thus, for a finite 
cluster of A^ sites one can compute the order-parameter Oi = Oi(K.,b,q) and the other 
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non-critical variable O2 = 02(K,b,q) (in the above examples O2 is the net magnetization 
m_(. + ■m_ of the antiferromagnet and the mean value of the square of the spin < Sf > in 
the Blume-Capel model), b is the field conjugated to the order-parameter Oi and q the 
respective field conjugated to the non-critical variable 02- Close to the transition just Oi is 
small so the expressions above are expanded for small b to give 

0, = F,{K,q)b, (58) 

O2 = O2{K,b = 0,q) = F2{K,q) . (59) 

Note that as b ^ one has Oi ^ and O2 7^ 0, irrespective of the value of g 7^ 0, whereas 
O2 — > when one further has q ^ 0. Similarly, for a smaller cluster of A^' sites one finds 

0[ = F[iK',q')b' , (60) 

O', = O',{K',b' = 0,q') = Fi{K',q'). (61) 

From now on we will omit the subscripts N and A^' and just denote the corresponding ther- 
modynamic functions for different clusters by unprimed and primed quantities, respectively. 
According to the ordinary MFRG one then finds 

Fi{K',q') = F^iK,q), (62) 

which, although independent of any exponent and rescaling factor, depends on the additional 
fields q and q'. It should be said that the above behavior is not quite general. Counter 
examples are, among others, disordered systems [0,0] and some quantum spin models ||5^ 



where, despite having a two parameter Hamiltonian, they present no non-critical variables 
being unnecessary the introduction of extra fields q and q'. As a result, Eq. (p2| ) is just 
dependent on K and K' and the procedure Ki —>■ K[ = R{Ki,r), for fixed values of the 
ratio r = K2/K1, outlined in the preamble of this subsection can easily be implemented. 
In studying the antiferromagnetic Ising model in the triangular lattice Slotte |S^ has 



proposed the mean field Ansatz where q and q' are self consistently obtained by requiring that 
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O2 = q' and O2 = q. Reasonable results, specially for the value of Kc have been achieved. 



This approach has also been used in the study of the Blume-Capel model [Q where results, 
not consistent with the expected ones (mainly by taking the smallest cluster sizes) have been 
obtained, including an unexpected phase transition at a finite critical temperature for the 
one-dimensional version (as seen in section II. A the exact result is expected to be obtained 
from smallest clusters at least for the one- dimensional model). 

A different choice, and equally natural, for the size dependence of the non-critical variable 
has been proposed by Plascak and Sa Barreto |^ and successfully applied to the Ashkin- 
Teller model. As this variable has its own size dependence which is not governed by finite 
size scaling at the transition it is assumed a renormalization group Ansatz in a way that 
q' = <1,0'2 = O2, or 

F^(K',g)=F2(K,g), (63) 

from which q = q' is obtained as a function of K and K' and inserted in Eq. ([6^). This 
renormalization group Ansatz has been proved to be better than the mean field one for 
the antiferromagnetic Ising model [Q and has been subsequently used in treating non- 



critical variables such as in the random Potts model |^. Concerning the Blume-Capel 
model defined in Eq. (p^), this improved formalism provides now exact results for its one- 
dimensional version and better values for the critical exponents and critical temperatures in 
dimensions d > 1 [^] when compared to those from the previous usual application of the 
method |58|. Besides, for the particular spin-1/2 case, the above formalism reproduces the 



early results of the Ising model since O2 =< Sf >= Oi = 1/4 is automatically satisfied. 

The renormalization group Ansatz allows one to obtain, in addition, an unambiguous 
estimate for the non-critical thermodynamic function O2 at criticality (note that the mean 
field Ansatz gives two different estimates O2 = q and O2 = q', respectively). However, it 
is not free from extracting results only in plausible invariant subsets given by fixed ratios 
r = K2/K1. So, complete analysis of the flux diagram is still not possible in such models. 
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2. two uncoupled order-parameters 

Some systems present two different and independent order-parameters which can be used 
to go beyond the hmitation (lacking of flux diagram) of the usual MFRG discussed so far. 
For example, in the symmetric Ashkin- Teller model |^3[ where two different spin variables 
exist on each lattice site, it is possible to identify two order-parameters in the problem. 
In a similar way, for the two-dimensional Ising model with crossing bonds one can select 
appropriate pairs of order-parameters depending on the region of the parameter space. In 
such models, for a finite block of A^ spins, one has Oi = Oi(K,6i) and O2 = 02(K,b2), 
where 61 and 62 are the symmetry breaking fields related to the order-parameters Oi and 
O2, respectively. Since for the uncoupled order-parameters Oi ^ when bi ^ (irrespective 
of the value of 62) and O2 —>■ when 62 ^ (irrespective of the value of 61) one has, to 
leading order in the b's, 

Oi = Fi(K)6i , O2 = F2{K)b2 . (64) 

By doing the same calculations for a smaller block of A^' sites and assuming that both 
order-parameters are described by the same critical exponent (which is usually the case) one 
gets 

F;(K')=Fi(K), (65) 

F^(K')=F2(K). (66) 



So, contrary to the procedure prescribed in Eqs. ([131 ) - ([18|) , in this case from equations (|65D 
and (^) one can obtain renormalization group flux diagrams of the usual type, furnishing 
isolated fixed points, critical exponents, universality classes, etc. Quite good results (includ- 
ing some exact values at some special points in the phase diagrams) have been obtained in 



the Ashkin- Teller model [^ and in the two-dimensional Ising model with crossing bonds 



through the above scheme. 
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3. two coupled order-parameters - EMFRG 



The independence of the order-parameters, as stated in Eq. (|MD, is not commonly verified 
and there are some models where they do not behave in such a simple form as discussed 
above. In the most general case one may have Oi = Oi(K, bi, 62) and O2 = C'2(K, 61, 62) so 
that, to leading order in 61 and 62, one gets 

Oi = Fn(K)6i + Fi2(K)62, (67) 

02 = F2i(K)6i + F22(K)62, (68) 

where both order-parameters depend on bi and 62- This general problem can be tackled by 
defining the vectors O = f^M and b = LM so that equations ( |57D and ( p^ ) can be written 
as 

O = Fb , (69) 

where the matrix F is given by 

(-^11 F12 \ 
F f) ■ '™' 

Similar expressions are obtained by considering a finite system A^' with less degrees of 
freedom. It is now easy to see that when the two order-parameters have the same critical 
exponents the corresponding vectors O and O' satisfy the scaling relation 

O' = L'^-y^O , (71) 

and the same for the symmetry breaking field vectors 

b' = L'^-y^h (72) 

if we take the usual MFRG. As a result, one ends up with the following matricial equation 

(F' - F)b' = . (73) 
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The above inatricial approach is a natural extension to the general problem of two-fold 
Hamiltonian parameters. This extended mean field renormalization group (EMFRG) repro- 
duces the previous cases as well as another particular approach proposed more recently by 



Likos and Maritan |^6[. Moreover, it can give the correct Bethe limit by using the smallest 
clusters in treating the mixed spin Ising model (an early MFRG study of the mixed spin- 
1/2-spin-l 1^ Ising model has failed to obtain the corresponding Bethe results due to a 



different mean field assumption in order to decouple the order parameters). Some procedures 
in treating of the matricial Eq. (|73D are now in order: 
(i) The simplest solution of Eq. ([73|) is given by 

det{F' - F) = . (74) 

From the above equation we do not have flow diagrams and the critical lines are obtained 
only for fixed values of r = K2/K1. This is a generalization of the corresponding one-order 
parameter where one has just F'(K') = -F(K). Now, by using the computed expressions 
of reference [^ for the two order-parameters of the mixed spin-l/2-spin-l Ising model one 
obtains the corresponding equation (173) which gives the same critical temperature from the 



Bethe approximation by using the clusters depicted in Fig. 1(a) and (b), as expected. 

(ii) If the two parameters are decoupled one has F^ = 0, i 7^ j for any cluster and 



another solution for Eq. (^S]) can be given by equating to zero all the matrix elements of 
F' - F, i.e., 

Fi, -Fn = and F^^ - F22 = , (75) 

which is the same result given by equations (|65D and (|6^). 

(iii) For the types of models considered by Likos and Maritan @,^] (in particular, 
models defined on a bipartite lattice with only two different ground states, a ferromagnetic 
and an antiferromagnetic one, separated by a stable borderline) one has Fu = F22 and 
F12 = F21 for any finite system. The last two relations reflect the additional symmetry 
Oi(K, bi, 62) = C'2(K, 62, ^1) exhibited by these order parameters. So, by equalling again to 
zero all the matrix elements of F' — F one gets 
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F[, -Fn = and F[^ - Fu = , (76) 



which are the same expressions as those proposed in reference [|6^ in a different context. 



One can thus see that the present matricial EMFRG is capable to generahze all the 
previous isolated treatments done on two-dimensional parameter Hamiltonians. One can 
also notice that: 

(a) in the general case where all the matrix elements of F are different from each other, 
Eq. (|73p can furnish four distinct recursion relations, i.e., at first sight more equations 



than parameters. However, it should be possible, in these cases, to augment the original 
Hamiltonian by a suitably chosen number of interactions until the number of couplings 
matches the number of recursion relations. From the flow diagram so obtained in this 
enlarged space one could get the corresponding flows in the original restricted domain by 
taking the appropriate subspace in the complete Hamiltonian space. 

(b) in addition, the present formalism can also be straightforwardly applied to models 
having more than two order-parameters simply by computing the dimensional enhanced 
vectors O and b and the corresponding matrix F for any finite cluster. 

Unfortunately, the extensions (a) and (b) above have not yet been tested in studying 
statistical mechanics systems from MFRG and one can not say, a priori, that this approach 
will succeed in giving the correct critical behavior in these cases, as is usual in ordinary real 
space renormalization group procedures. 

D. Dynamic MFRG (DMFRG) 

So far, just equilibrium static properties have been treated through MFRG. In this 
subsection we discuss the extension of the mean field renormalization group ideas to non- 
equilibrium phenomena. Despite the simplicity of the formalism, the MFRG has been 
applied to a very few number of systems out of the thermodynamic equilibrium. It was 



employed for the first time by Indekeu, Stella and Zhang p3 to study the dynamics of the 



kinetic Ising model with single spin-flip Glauber transitions |T^ near the equilibrium states. 
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It has also been used to treat the dynamic critical properties of quantum spin systems [|7T . 
As the extensions to treat the dynamics of quantum models is (at least in principle) straight- 
forward to be done, we present below just the formalism for treating kinetic classical models 
with special emphasis to the Ising model by considering larger clusters than in the original 



work of Indekeu et al [pU|. Extensions to other problems are rather easy to be performed 
and will be summarised in section V. 

In the dynamic approach of the MFRG one starts from the scaling relation of the mag- 
netization, close to equilibrium and criticality. For sufficiently long times t we expect the 
magnetization of the infinite system to scale as ||72| 

m(e, H, t) = r'^+y^mif/'^e, i^^H, rH), (77) 

which is a generalization of Eq. (|^) where z is the corresponding dynamic critical exponent. 
Let S = {Si, S2, . ■ . , Si, . . . , Sn), with Sj = ±1, represent a state of the finite system 
with N spins, and P{S,t) the probability of finding the system in the state S at instant t. 
The time evolution of P{S,t) is given by the master equation 

^^^ = J:[PiS\t)W,{S')-P{S,t)W.{S)] , (78) 

where Wi{S) is the transition rate, per unit time, to fiip the spin Si, and 5* = 
(^i, 5*2, ... , —Si, . . . , Sn)- The first sum in the above equation takes into account all transi- 
tions to the state S and the second sum all the transitions out from the state S. If f{S) is 
a given function of state S, we can evaluate its average value by 

<f{S)>=j:fiS')P{S',t), (79) 

S' 

and its time evolution through the expression 

By taking the function of state as being the local magnetization of the spin i, we can easily 
write 
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'^^^ = -2 < S,W,{S) > . (81) 

Let m]\f{t) and mi\i/{t') be the time dependent magnetizations of two finite clusters with 

N and N' spins, respectively, i.e., 

1 ^ 
mN{t) = j;Tj:<Si>, (82) 

1 ^' 



where the time dependence comes from < Si > through Eq. (pTI). Following the original 
idea of the MFRG , each spin Si of the border of the cluster couples with an infinitesimal 
symmetry breaking field b(t). The equations of motion for mN{t) and m]^i{t') are computed 
from Eqs. (|T|)-(|81]) by using |7| 

Wi{S) = \[^-Si tanh(ir ^ 5, + aKhN)] , (84) 

where q is the number of first neighbors of Si, which are external to the cluster, and 

Yl qijkit)SiSjSk + . . .) , (85) 

i,j,k 

where mi(t) =< Si >, Pij(t) =< SiSj >, qijkif) =< SiSjSk >, etc. The transition rate Wi{S) 
is chosen to satisfy the detailed balance condition P{S\ t)Wi{S^) = P{S, t)Wi{S) in the limit 
oi t —>■ oo. This is a sufficient condition to the equilibrium be attained, i.e., dP{S, t)/dt = 
as t — > oo with the expected probability distribution Piy{S,t -^ oo) oc e~^^^'^\ where E{S) 
is the energy of the state S. The expression for the probability P/v(5,t), Eq. (p5|) , is an 
exact one and can be worked out just for rather small systems. 

In general, for small clusters a.t H = and linearized in b]\}{t) and hj^i{t') the time 
derivatives of the order parameters can be written as 
dniNit) 



dt 
dniN'it') 



AN{K)mN{t) + BN{K)bN{t) , (86) 

ANiK')mN,{t') + BN\K')hN\t') , (87) 



dt' 
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where the expressions for An{K), An'{K'), Bm{K) and Bn'{K') are computed from Eqs. 

m-m- 



According to the strategy of MFRG, we impose now the following scaling relations for 
the magnetizations and symmetry breaking fields: 

mN>{K',Q,rH-hN>{rH)) =i'^-y"mN{K,0,t;bN{t)) , (88) 

and 

bj^,{rH) = i''-y"bN{t) , (89) 

where i = i§)^ and t' = i-H with A^' < A^. 

Taking the derivative of Eq. ( pS]) with respect to t, and using Eqs. (p6D, {^7\) and (p9D, 



we arrive at the following recursion relations for the A's and B's coefficients: 

Ar,iK')=tAr,{K) , (90) 

B^iK')=eBN{K) . (91) 

The solution of the above system of equations provides the non-trivial fixed point {K = 
K' = Kc) and furnishes also the value of the dynamic critical exponent z. 

As a simple example, let us take the one and two spin clusters depicted in Fig. 1(a) an 
(b) for this dynamic treatment in the ci-dimensional Ising model. We have then 

^^ = -2 E ^4[1 - S,t^nh{cK'b,)]^[l + m^S^] , (92) 

{Si} 

^^ = -2 Y. S,l[l-S,te.nh{KS2 + {c-l)Kb2)]x 

x-[l + m2(t)5i + m2(t)52+pi2(t)5i52] , (93) 

where in the latter relation we have used the fact that both spins in the two site cluster 
of Fig. 1(b) have the same time dependence. Performing the sum over the corresponding 
states in the two relations above one gets 
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dt 
from which one has 



^^ = -m,{t) + cK'hit) , (94) 

{l-tanhK)m2{t) + {c-l)K{l-tanh'^K)b2{t) , (95) 



1 =r(l-tanhi^) , (96) 

cK' = t{c - l)K{l - tanh^ K) . (97) 



By taking the ratio of Eq. (^) to Eq. (|96|) it is easy to see that the non-trivial fixed 



point and the related thermal critical exponent are the same as those obtained in the static 
procedure of section IIA (as well as the magnetic critical exponent yu if one includes the 
external field if), as expected. Estimates of the corresponding dynamic critical exponent z 
are also readily achieved from above equations (see Table 0). 

In Table we show the results obtained in the study of the Ising model for clusters of 
up to 9 spins in two dimensions (a cluster with non-equivalent boundary sites which has not 
been previously considered in the dynamic treatment of the model) and with up to 8 spins 
in three dimensions. One can see from this Table that the results are consistent with those 
from other approaches. However, some difficulties arise in considering bigger clusters in the 
present case: 

(a) when we consider a cluster with four spins, by using Eq. (^) it appears correlations 
between pairs of spins which are located at the diagonals of the cluster; 

(b) we must also be careful in the calculation of the mean values < Sj >, because for 
clusters with more than four spins they are not all equivalent. It should be stressed that 
this non-equivalence among spins inside the finite cluster for A^ > 4 in two- dimensions (and 
A^ > 8 in three-dimensions) implies in the necessity to diagonalize a system of equations in 
order to obtain an expression like that of Eq.(^), since in this case vn^qif) is no longer an 
eigenmode of the dynamics but is coupled to higher spin correlations; 

(c) besides, another great difficulty arises when we try to increase the size of the clusters 
in this non-equilibrium approach of the MFRG. Although the method does not produce 
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couplings of longer range in space, it is necessary to take account of the higher than two 
spin correlations, which naturally appear for bigger clusters. Then, we need to use methods 
of equilibrium statistical mechanics to evaluate these higher order correlations, although we 
are interested only in the determination of the equation of motion for the mean magnetization 
of the cluster. 

As a result, the dynamic MFRG must be applied only to small clusters. 

III. EXACT RESULTS FROM MFRG 

Besides the great qualitatively success achieved by applying the MFRG to a great variety 
of statistical systems (even by taking the simplest choice for the clusters) some expected 
rigorous quantities are obtained in some particular cases, as we have already seen in the 
previous section. That is what happens, for instance, when treating the one-dimensional 
spin-1/2 Ising model |2^, the two-dimensional spin-1/2 isotropic Heisenberg model |]75[ and 



the spin-S one-dimensional Blume-Capel model |^2|. While the exact critical temperature 
[Tc = 0) is reproduced in all models above for clusters of any size, in some of them the exact 
critical exponent is only obtained in the limit of infinite cluster (except for the magnetic 
exponent i/h in the Ising model and the thermal exponent z/ in the d = 2 Heisenberg model, 
where they are exact even for the smallest clusters). 

One also obtains expected quantitative values at some special points in the global phase 
diagram of the Ashkin- Teller model [^,0 and the Ising model with crossing bonds |]65| , |66| . 
Exact results are also achieved by applying the method in the study of the geometrical 
problem of directed self-avoiding random walk in two dimensions and the Ising model on 



the Bethe lattice [^]. These rigorous results reproduced by the MFRG in so many different 
systems such as classical, quantum and geometrical statistical problems assign indeed to the 
method a high reliability when studying more complex models. 
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IV. STATIC PROBLEMS APPLICATIONS 



The previous two sections give us an idea of the performance and potentiahty of the 
phenomenological mean field renormahzation group when used to treat some particular 
models. In this and in the next section we will present, in a summarized way, the applications 
of the method (in its various forms) to several different statistical mechanical systems. It is 
not the purpose here to discuss the results obtained for every studied system (which amounts 
to over 120 papers). Rather, just an update of the treated models by employing the method 
will be presented. Moreover, this section is intended to be self-contained in the sense that 
all the references already commented in the text is referenced here again when listing the 
corresponding model where it has been used. 

The Ising model has been extensively studied and has also been a test system for all 
proposed procedures involving the MFRG approach [^ , ^ , ^9| - |5Tl , |76| - f7^] , including a com- 



bination with the decimation procedure ||79|,|80[|, the effects of reaction |81] and symmetry 



dependent fields [g^], and treatments of some extended models such as: antiferromagnetic 
model on triangular |^ and square lattices with external field [0]; anisotropic Ising model 
on a triangular lattice [|68|; surface effects in pure and random semi-infinite models ||83|-|93 



the Ising system on a compressible lattice (Domb model) ||9^; and the model with crossing 

bonds in two- dimensions [^,^,^. Additional diluted versions of the Ising model have been 

treated by considering the effects of different kinds of random dilution including spin-glass 

55| , |8llj96Hl09|] , dilution in the antiferromagnetic model in an external field ||110|| and random 



field systems [p5| , |lll 



Geometrical critical phenomena and percolation have been studied by De'Bell |^,@] and 

das Neves and Kamphorst Leal da Silva ||112|| , and some exact results have been achieved 

48| . The method has also been used in the study of directed percolation hypothesis for 



stochastic cellular automata |113|, lattice gas model |114|, and nematic ordered states at 



low temperatures |p.l5| , [116 



Other related systems include: Ashkin- Teller model p9| , |6^ , pTr| ; coupled spin-1/2 Ising 
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models ||118|| ; ANNNI model ||119|| ; mixed spin-1/2 spin-1 model |]67| , |120|| ; different versions 
of pure and random Potts models on infinite and semi- infinite lattices [ |61| , |121| - [127|| ; the Z4 
spin model ||128|| ; Hamiltonian version of the two-dimensional Z{q) symmetric spin models 
|129| , |130|| ; quenched and annealed random-bond D-vector models ||131|| ; the planar random 
anisotropic model [ p.32|| ; g-state clock spin-glass models | |133| | ; classical pB| , |52| , pr51| - |136| ] and 
quantum random versions of the Blume-Capel model [|137|| ; the Blume- Emmery-Griffiths 



model ||138|| ; the finite temperature SU(2) lattice gauge theory at the strong coupling limit 
139|| ; and classical XY and Heisenberg models ||140| . 



A great variety of quantum spin systems have already been studied from the MFRG 
approach. Among them we have: anisotropic Heisenberg model |[75|J141| - [144|| ; anisotropic 
Heisenberg model in a transverse field ||145|| ; pure and random spin-1/2 [^ , |146| , p!47[l and 
spin>l/2 ||148| - |15l| transverse Ising models; quantum Ising model with annealed antiferro- 
magnetic bond randomness ||152|| ; Ising model in random transverse field [ p.5iyi53H1571 



ran- 



dom mixture of Ising and Heisenberg models |[158|J159|| ; spin-1 anisotropic Heisenberg chain 
16U| and spin-1 Ising model in a transverse field ||161| ; quantum models applied to gran- 
ular superconductors ||162|| ; and anisotropic Heisenberg model with Dzyaloshinskii-Moryia 
interaction ||163[ . 

More recent extensions of the method include the study of the MFRG approach with the 
Tsallis Statistics ||164|| and non ambiguous location of multicritical points ||136| . 



V. DYNAMIC PROBLEMS APPLICATIONS 



Although fewer in number, the DMFRG has been applied in the study of dynamics and 
non-equilibrium properties of some interesting systems. For instance, the classical approach 
of Indekeu et al [^ has been extended in the treatment of the intrinsic dynamic critical 
properties of quantum spin models [|7T|. In this case, good results for the critical exponent 
z have been obtained when compared to others from more laborious approaches. 

It was also applied to non-equilibrium problems arising from the competition between 
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microscopic mechanisms. This is the case, for instance, of the Ising model subject to two 
locally competing temperatures ||165| , |166|| . In this model, for each temperature, the corre- 



sponding transition rate satisfies the principle of detailed balance. When the two processes 
are considered at the same time, a continuous phase transition is observed between steady 
states. By using clusters of one and two spins Marques ||165| | was able to determine the 



phase diagram of the model in the temperature versus gradient of temperature plane. She 
also obtained values for the critical exponent v of this non-equilibrium model in two and 
three dimensions, and they compare very well with the values of the corresponding equi- 
librium Ising model. On the other hand, if one of the temperatures becomes negative, the 



ferromagnetic system can display an antiferromagnetic order [p.66|| . The heat bath at neg- 
ative temperature simulates a flux of energy into the system. The MFRG was also applied 
to other non-equilibrium models which exhibit steady states: non-equilibrium Ising model 



with competing dynamics ||167|| , critical behavior of non-equilibrium 3-state systems ||168 



stochastic cellular automata in one and two dimensions [ |113|| , competition between diffusion 
and anihilation of many particles ||169|| , dilution in the contact process | |170|| , critical surface 



behavior of models with one absorbing state | p.71|| and one-dimensional models with multiple 
absorbing configurations | |172 . 



This methodology has also been used to study the dynamic critical behavior of a semi- 
infinite system of spins in a simple cubic lattice with nearest-neighbor ferromagnetic in- 
teractions. As this approach has not been published in the literature and presents some 
similarities with the one discussed in section II.C.3 we outline below the general procedure. 
It is assumed that in the surface plane of the semi-infinite system the exchange coupling is 
given by Ks-, while all the other couplings are given by K. For p = Ks/K greater than a 
value pc the surface orders before the bulk while for p < pc the surface orders when the bulk 
does. As an example, in the one-site cluster of the MFRG, we take one spin in the i — th 
plane, and it interacts with the symmetry breaking fields: b\{t) of the i — th plane, and 
b]^^(t) and b\^^(t) of the adjacent planes. Here, and from now on, the superscript designates 
the particular plane away from the i = 1 surface. As is usual, these fields are assumed to 
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be very small, that is, we are considering the problem in the neighborhood of the surface 
phase transition. The equation of motion for the magnetization taking the one-site cluster, 
in each plane, is then given by 

^^ = -m\ it) + AK'sh\ it) + K'h\ it) , (98) 

for the magnetization on the surface i = 1 and, for i > 2, 

^"^- ''^ ^ _^^ (^) ^ K'[b\-\t) + 4b\{t) + b\+\t)] . (99) 



dt 
Eq. (|99|) is a generalization of Eq. (|94D for the bulk magnetization with different sym 



metry breaking fields in different planes and Eq. (|98D is a particular case of relation (|99D 
for the surface magnetization with intra-plane interactions Ks- Similarly, the generalized 
expressions for the two spin cluster contained in each plane read 

dml{t) 



dt 



-(1 - t&nh Ks)ml{t) + (1 - tanh^ Kg) [3 Ksbl{t) + Kbl{t)] , (100) 



^^^^1^ = -(1 - tSinhK)mi{t) + (1 - tanh^ K)[KUi^ + ?,KU^{t) + KU+\t)\ . (101) 

In general, the magnetization for the planes i taking a small block of N spins inside each 
layer can be written as 

dm^ (f) ^ 

-^ = -Al,{KsWN{t) + Y.B'AK,Ks)Kit) , (102) 

dm^ (f) *+^ 

^-^ = -AV(ir)mV(t) + Y: B^AKWNit) , ^>2, (103) 

where the explicit dependence of the coefficients AJ^ and B^j on K and Ks has been given. 
By considering now two clusters of A^ and N' spins and the scaling relations given by 
equations ( |88| ) and ( |89| ) for the magnetizations on each plane (and introducing adequately 
the exponents i/hs and zs for the surface magnetization and i/h and z for the bulk quantities) 
one finds, after taking the derivative with respect to the time t, a set of linear equations 
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relating the magnetizations and the corresponding symmetry breaking fields of different 
layers. The set of linear equations for the magnetizations can be written as 

AnMn = An'Mn' , (104) 

where A^ at is a column vector composed by the layer magnetizations m]^, mj^, m\, ..., and 
An is a diagonal matrix whose elements are straight related to A^. The same holds for the 
smaller system. A solution of (|104|) is obtained by equating all the terms of the diagonal of 



the matrices and, as all terms for i > 2 are the same, results basically in 

A\{Ks) = r^'A],,{K's) , Aiv(i^) = i-'A'AK') , (105) 

which gives a relation among Ks, K, z and zs at criticality. The set of linear equations for 
the remaining symmetry breaking fields can be put in the following form 

Bh = , (106) 

where b is the column vector composed by the symmetry breaking fields b\f, bj^, 6^, ..., 
and i3 is a three-diagonal matrix whose elements, using relations ( |1U5D , are just functions of 
Ks and K at criticality. It is easy to see that the critical coupling to give an ordered surface 
over a paramagnetic bulk phase is found when the determinant of matrix B vanishes. One 
then obtains an equation for the surface critical temperature as a function of p = Ks/K, 
as is usual in the MFRG approach of two parameter systems. The previous bulk dynamic 
properties are reproduced, as expected, for i :» 1. 



For the simplest finite systems of Fig. 1 it is not difficult to show from Eqs. (p8| ) - ([T0lD 
and (|T0|)-(|106D that 

tanh Ksc .,„„s 

1 - ManhKs. = ^ ' ^'^^^ 

for the surface critical temperature Ksc as a function of the ratio p and 

1 - tanhKs = l'^' , 1 - tanhK = l'^ , (108) 
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from which the dynamic critical exponents are obtained. The numerical results from above 
equations are shown in Table [V^ together with those obtained by taking the plaquette of 
four spins. We can note that the DMFRG is easily extended to study the critical dynamic 
phenomena on surfaces and, in particular, furnishes a value for pc = 1.48 comparable to 
that from Monte Carlo simulations pc = 1.52 | |1 73 1 , while the mean field result is pc = 1.25 



36| , |174|| . It should also be noticed that in the present approach with rather small clusters the 



surface magnetic exponent yns is indeterminate since it cancels out in the renormalization 
group equations. 

VI. RELATED PHENOMENOLOGICAL RG 

Different phenomenological renormalization group procedures based on exact calcula- 
tions in finite systems have also been proposed in the literature. In this section some of 
them will be discussed, namely, the effective field renormalization and finite size scaling 
renormalization groups. Few other methods, which are less used or applied in approaches 
other than the renormalization group scheme, will be referenced in the next section. Here, 
as in the MFRG case, we are mainly concerned in describing the methods themselves and 
their interrelations and similarities, and just a list of their applications in different statistical 
models will be presented. 

A. Effective Field Renormalization Group (EFRG) 

The MFRG in its several formulations computes exactly the magnetization (order pa- 
rameter) for each finite cluster of spins according to the canonical distribution (plQl). An 
alternative way of obtaining this order parameter for a Hamiltonian model system 7i has 
been proposed by some authors ||175| - |177|| by employing an effective field theory based on 



the exact generalized Callen-Suzuki identity ||178yi79| 



<«")-{^^ai^>' <-' 
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where the partial trace Tr is taken over the set of n spins variables specified by a finite 
system Hamiltonian 7i„, 0„ is the corresponding order parameter (or, in general, any other 
function of all the n spins of the cluster) and < • • ■ > indicates the usual canonical thermal 
average taken over the ensemble defined by the complete Hamiltonian 7i. The idea is simply 
to replace ?^ by a cluster of n = A^ spins surrounded by fixed magnetizations at values b. 
In this way the order parameter m = O in Eq. ( |109| ) can be computed by employing the 
exponential operator technique [|180|| resulting, for 6 << 1, in an equation of a similar form 



as before for H = 0, namely, 

m^ = /^(K)6^ . (110) 

Doing the same for a smaller system (A^' < A^) and using the mean field renormalization 
group assumptions given by equations (pO|) and (^) one obtains the EFRG recursion relation 

/^,(K') = /Ar(K), (111) 

from which phase diagrams and estimates of critical exponents are computed. 

In order to illustrate the method let us consider the spin-1/2 Ising model and n = 1 in 
Eq. (11091) . The Callen-Suzuki identity ([109D then reads 

mi=<Si>=/tanh{K'Y,Sj)\ , (112) 

where the sum is over all nearest neighbors spins of Si. By using the exponential operator 
technique ||180 



e'""- fix) = fix + 5), ^- = ^' (113) 



the Eq. (|112|) can be written as 



mi = (e'^^- tanh(x) U=o) , 6 = K'J2Sj , (114) 

and, as the hyperbolic tangent does not depend on any spin configuration one has 
mi = lfl e^'^^^- \ tanh(a;) |^=o= ( H [cosh(X'L)^) + Sj sinhiK'D^)] \ tanhx |^=o , (115) 
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where in the last expression we have used the van der Waerden identity for the two state 
system 



^aSi 



cosh(a) + S'iSinh(a), ^j = ±1 . (116) 



The above equations are still exact relations and computation of irti will commonly require 
some approximations such as decoupling the mean values of product of spins in product of 



spin mean values ||180|| . We can also see from Eq. ( p,15| ) that the mean value mi is strongly 
dependent on the number of first neighbors c. By taking the one spin cluster A^' = 1 of Fig. 
1(a) we see that Sj = hi for all j. The order-parameter for the one-dimensional lattice Eq. 
assumes then the form 

mi = ([cosh(ir'L)^.) + 6isinh(K'L)^)]2\tanh(a;) U=o • (117) 



Since sinh(ir'L)^.) cosh(ir'L)^.) tanh(x)|^.=o = |tanh(2ir'), sinh2(ir'D^) tanh(x)|^.=o = 
and cosh (K'D^) tanh(a;)|^=o = one has for a homogeneous system 

mi = tanh(2is:') bi . (118) 

It is interesting to note that the usual mean field approach mi = bi gives the exact result 
Tc = even for the one spin cluster. In fact, this approach reproduces all the expected 
results for the thermodynamic properties, such as magnetic susceptibility and specific heat, 
of the linear chain in zero external field [|181|| and reflects the fact that, besides appearing 



at most just pair correlation functions in the thermodynamic functions, the auto-correlation 
< Sf >= 1 has been taken into account exactly in the van der Waerden relation ([1161) . 
Analogously, for the two-dimensional square lattice one gets [[180| 



mi = i[tanh(4K) + 2tanh(2K)]6i , (119) 

which furnishes the mean field critical temperature Kc = 0.324, and for the three-dimensional 
cubic lattice 

3 

mi = — [tanh(6K') + 4tanh(4ir') + 5tanh(2K')] bi, (120) 
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giving Kc = 0.197 [ p.80|| . As usual, from the above simple approach one always obtains 



classical critical exponents (as can be seen by expanding further up to the order bf] 
For a two spin cluster n = 2 Eq. ( p.09| ) reduces to ||182|] 



,„, 4 (5. . 5.) ^ (fie-.- fie-.-') „.i,(.^:;";;^:i„._,) . ^m 

where Dy = 4- and the sum in j is over the first neighbors of the spin 1, and the sum in j' 
is over the first neighbors of the spin 2. It is now easy, though rather lengthy, to compute 
the expressions for m2 taking the cluster of Fig. 1(b) on different lattices. We quote below 
just the results for the one- and two(square)-dimensional lattices which are, respectively, 

_ 3 / /(6) 4/(4) 2/(2) 3/(2) \ 

""' " 24 1^^(6) + e-2^ + ^(4) + e-^^g{2) + g{2) + e~^^ g{A) + ^(2) + e-^^ } '^ ^ > 

where f{u) = sinhwi^ and g{u) = coshuK. A complete set of coefficients for other types 
of lattices taking the smallest clusters can be found in reference [ p. 77|| . Eq. ( |122| ) gives also 



Tc = for the one- dimensional case and one gets K^ = 0.331 and K^ = 0.198 for d = 2 and 
d = 3, respectively 



The EFRG obtained from Eqs. ([TToD , ([TTll) , ( ITTSl ) and (|T2|) reproduce the exact one- 



dimensional results, including the expected thermal critical exponent u even for the smallest 
systems. In general, the corresponding renormalization recursion relations for dimensions 



d > 1 give better values than the MFRG, as can be seen in Table [VII] , where a comparison 
is made on square and simple cubic Ising models. In addition, the two spin block is, within 
this formalism, able to distinguish the two dimensional triangular lattice from the simple 
cubic one (both with the same coordination number c = 6) regarding the critical temperature 



{Kc = 0.244 for the triangular lattice and the exact value is Kc = 0.275 |^). This distinction 
comes from the fact that on a triangular lattice the total number of surrounding mean field 
sites for N = 2 cluster is 8 and in the simple cubic one is 10. While Eq. ( |121[ ) depends on this 
value, in the MFRG approach expressed by Eq. (|38D it only matters the number of neighbors 
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on each site and, in some sense, the common neighbors on the triangular lattice is counted 
twice. The price paid, however, for such improvements in rather small systems taking the 
EFRG procedure clearly reflects the more elaborate calculations needed to obtaining the 
order parameters. As a result, contrary to the MFRG, this drastically limits the size of the 
employed clusters in such a way that the bigger system used in literature is A^ = 4. 

Besides pure and diluted Ising models ||175H177| , |183H186|| the EFRG has also been apphed 



in the study of the Ashkin- Teller model [p.76|] , surface criticality in semi-inflnite Ising fer 



romagnets ||187|| , classical XY and Heisenberg models [ |14CI|] , 0{n) vector ferromagnetic and 



antiferromagnetic models [|188|| , the transverse Ising model [|189| , |19(]|| , and quantum spin-1/2 



anisotropic Heisenberg models ||191| , |192|| . The same qualitative results have been achieved 
in all applications. 

B. Finite Size Scaling Renormalization Group (FSSRG) 

The finite size scaling renormalization group (FSSRG) has been originally proposed by 
de Oliveira p5|,pP| to treat Hamiltonian systems having two energy terms. The main idea 



of the method is to consider quantities having zero anomalous dimension, i.e., = in 
equation (0). For the simple Ising model in an external magnetic field defined in Eq. (|l]) 
such quantities can be given by 

Q{K,H) = /s^gn(^Y.S^]) , (124) 



and 

R{K, H) = l^sign i^f. ^A sign i^f. ^A ^ , (125) 

where N = L'^ is the number of spins on a d-dimensional lattice and sign{x) = —1,0,1, 
whether x < 0, x = 0, x > 0, respectively. Si and S2 are two parallel surfaces with Ns^ = Ng^ 
and separated by a distance L/2 apart when taking periodic boundary conditions. For finite 
lattices with open boundary conditions the surfaces Si and S2 can be taken as being the top 
and the bottom hypersurfaces of the corresponding hypercube, respectively. 
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Let us consider first the quantity Q and see that it has zero anomalous dimension for an 
infinite system. While the magnetization m = (ij^Y^fSAj is zero for T > Tc and behaves 
as 

m^e'^ , (126) 

when H —^ 0^ and T — > T~, the quantity Q = +1 for T < T^ and is zero above T^ meaning 
that Q = e° and, from Eq. (0), = 0. Regarding now the second quantity R we have 
that the probability of finding -^J^i "^ Si and j^J^i "^ Si with the same sign is greater 
than finding them with opposite signs since, according to the weak version of the Griffiths 
inequahty [ |193|| , one has (SiSj) > for any pair of spins on the lattice. Thus, for the same 
reasoning, we also have a zero anomalous exponent for R. 

By taking two finite systems of sizes N, N' and computing the above quantities for both 
clusters one then gets the renormalization recursion relations 

QN>iK',H') = QN{K,H) , (127) 

and 

Rn'{K', H') = Rn{K, H) . (128) 

From equations above the complete renormalization flow in the K — H plane can be 
exploited. It turns out that for the Ising model the FSSRG gives the quantitative correct 
behavior as A^ > A^' ^ oo for all the fixed points {K -^ cxd, H = 0), {K = 0,H = 0) and the 
Ising critical point {Kc,H = 0), while expected qualitative behavior is achieved for finite 
systems. For instance, taking A^ = 4 and A^' = 2, which are the smallest possible systems 
allowed by this approach, one readily gets 

smh{2H') _ e^^ smh{AH) + Asmh{2H) 



cosh(2i/') - e-2^' _ e^-^cosh(4A') - 1 

7^ — 7^ ) (130) 
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where 



Zo = cosh2/f' + e^2-^" 



(131) 



Z4 = 2 + e"^^ + 4 cosh(2i7) + e^^ cosh(4/7) , (132) 

from which one obtains the values quoted in Table |Vlll| for the Ising critical point. Results 
from analytical calculations on bigger lattices are also listed in that table. 

However, contrary to the previous approaches, the FSSRG method has the great advan- 
tage of allowing Monte Carlo simulations for bigger finite systems in order to obtain the 
desired quantities Q and R . Such Monte Carlo simulations are not so easily implemented 
and has not yet been done in the mean field like renormalization group treatments discussed 
herein (apart from a theory of mean field Monte Carlo simulation proposed by Netz and 
Berker in a different context ||194j - |196|] ). These results are also shown in Table |VTT| and, even 
taking into account the rather small computational effort, they are quite close to the exact 
ones. The FSSRG can also be used to treat other systems with different kinds of energy and 
coupling constants by conveniently choosing the quantities Q and R and even, in principle, 
be extended to study more than two parameter Hamiltonians, such as flux diagrams in three 
dimensions. 

Apart from the pure Ising model [p8ip9ip!97|l , the FSSRG has also been applied to the 



diluted Ising model |Q, the spin-1 and spin-3/2 Blume-Capel model ||198|| , and the q-state 



Potts model [[L99|| . An extension to the dynamic critical behavior (for short and long times) 
of the Ising P00|| and Potts model |P01|| has also been done through Monte Carlo simulations 
on finite systems. In this case, quite good accurate results for the dynamic critical exponents 
z have been achieved. 

It should also be mentioned that the FSSRG approach is identical to the phenomeno- 
logical renormalization group proposed by Nightingale |]19l when in the former one con- 



siders clusters consisting of infinite strips with finite width ||202|| . In order to see this 
equivalence let us take a strip with finite width L and length n^ -^ 00 dX zero exter- 
nal field (for simplicity). Correspondingly, the transfer matrix will be of order 2^ x 2^. 
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In this case the quantity Ql-, computed from the transfer matrix, is zero and one gets 
Qli{K\ H' = 0) = Ql{K, H = 0) = ioi any values of L and L' < L. On the other hand, 
the quantity Rl is just a correlation function between the states of two columns separated 



by a distance ni/2 apart, which can be written as ||203|| 



,<\"i/2 



Rl{K) = A, l^^j , (133) 

where Al is a function of K but independent of ul and A^ and A^ are the greatest and 
second greatest eigenvalues of the corresponding transfer matrix, respectively. For two strips 
with L' < L one then has 






^^'(t# =^Mt# ' (134) 



where ul, uli -^ oo. Taking now the logarithm of both sides of the above equation, dividing 
by ny jl and recalling that the scaling factor I = L/L' can also be obtained from 

^ ' ^ * ' ' (135) 



(from which we have ^ = -^) we finally get 



ih-a: 



L'ln -3^ =Lln 1^ , (136) 



which is the same as Eq. (p!2| ) since ^^^ = In ( ^ j (the same for L'). 

It is interesting that two completely different phenomenological renormalization group 
procedures have the same (quite good) results for the properties of one parameter criti- 
cal models. One can think of the FSSRG as a generalization of the original Nightingale 
renormalization group for completely finite systems (a fact that, in some way, explains the 
accuracy of the FSSRG critical values, specially in the Ising case). On the other hand, flux 
diagrams for two-dimensional parameter Hamiltonians (for instance, the Blume-Capel model 



| 204|| ) taking infinite strips within FSSRG can be readily obtained (which is not possible 



from the Nightingale procedure). 
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Another real space renormalization group approach preserving the two spin correlation 
function in finite clusters proposed by Tsallis et al [0,^] has also been shown to be con- 
nected to the Nightingale procedure. However, it seems rather difficult to implement the 
calculations by taking infinite self-dual clusters in the former method. 

VII. FINAL REMARKS 

In the above discussion we have been mainly concerned with renormalization group pro- 
cedures based on finite size scaling hypothesis having a relatively wide application in rather 
different statistical mechanics problems. There are, however, some other renormalization 
group approaches based on ffnite lattice calculations. One of them, still close in spirit to the 



ones we have been described so far, is the new mean ffeld renormalization group P05|| trans- 
formation (NMFRG). To brieffy illustrate this method let us recall Eq. (^) in the mean 
ffeld approximation ttin = b. Close to the mean ffeld transition T^ the order parameter can 
be expanded as 

m^ = Fn{K)4 + GN{K)He~/ , (137) 

where e = {K — K^)/K, K^ is the mean ffeld critical temperature and (3* and 7* are the 
usual mean ffeld exponents {(3* = 1/2 and 7* = 1). At K = K^ the functions Fn{K) 
and Gn{K) are called anomaly coefficients, introduced by Suzuki ||206| , |207|| and from which 



non-classical exponents are also achieved (the so-called coherent anomaly method - CAM). 
In the present renormalization group context one expands the order parameter for two 
different clusters N' < N and assumes a scaling relation of the form given in Eq. ( pO]) for 
the approximate order parameters and the same relation for the quantities e^ and e^/, i.e., 
ejY' = i'^~'^"e^ . So, in the same way as previously done for the MFRG one then gets 

Fn>{K') = Fn{K) , (138) 

from which the critical ffxed point Kc = K' = K is obtained, as well as the thermal exponent 
z/, whereas the corresponding magnetic exponent is given by 
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<.«.(£^)-^-^. (139, 

It is apparent that although the anomaly coefficients are based on the classical values (3* 
and 7*, non-classical exponents are extracted from the linearization of the recursion relation 



( p.38|) around the fixed point. This approach gives quite good results for the Ising model 
as well as for geometrical phase transitions [p05|| , and can also be extended by using the 



effective field theory through the use of Callen-Suzuki identity ||208|| . It is also verified that 



the larger the value of N and A^', the better the approximation since more fluctuations are 
included in the mean field calculations within each cluster. 

Just to mention, another approach, used in a different context, is the density ma- 
trix renormalization group p09| , pi(]| ] which gives accurate numerical results in studying 
one-dimensional quantum lattice models [|211|] , the two-dimensional Ising model [|212|| and 



strongly correlated electron systems [pl3|| . Nevertheless, it is clear that other phenomeno- 
logical renormalization group proposals can still be done by exploiting relation ([T0|) for 
quantities where is known in the infinite system. 

Concerning the present schemes there are still some open problems whose study would 
be natural extensions of the methods discussed herein. Among them, for instance, we have: 

(i) there is no unified approach (as in the MFRG case) to the SBMFRG (and its cor- 
ner field partner) in order to get complete flux diagrams for models with more than one 
Hamiltonian parameter space. In fact, just for the Ashkin- Teller model a flow diagram in 
the Hamiltonian parameter space has been obtained through the SBMFRG ||64|| ; 



(ii) there is also a lack of a SBMFRG scheme to treat critical dynamics by considering 
three clusters. Only the two cluster approach has been implemented in the MFRG to 
dynamics; 

(iii) it would be interesting to extend the EFRG procedure to account for obtaining flow 
diagrams; 

(iv) flow diagrams in more than two Hamiltonian parameter space has not yet been 
considered by using neither the MFRG nor the FSSRG approach; 
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(v) finally, Table |^ shows roughly a picture of the models treated by the MFRG approach 
presented in sections III-V. It is also presented the applications of the other procedures 
discussed herein. So, besides new systems to be treated and studied from the present 
approaches, all the lacunae in Table ^ are straight extensions to be done on this subject. 
Perhaps, some readers will find it easy to readily fill in some gaps. 
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TABLES 

TABLE I. 
Critical values for the d = 2 Ising model on a square lattice according to the ordinary MFRG 
(two clusters) and SBMFRG (three clusters). The two entries for the critical indices in the 
SBMFRG correspond to the values given from equations (^) and (^Tj), respectively '^. 



MFRG 


N,N' 


Kc 


jy 


VH 


VHS 


^VH 


4,1 


0.361 


1.45 


1.50 




2.17 


9,4 


0.381 


1.28 


1.57 




2.01 


16, 9 


0.393 


1.22 


1.60 




1.95 


25,16 


0.401 


1.19 


1.62 




1.93 


SBMFRG 


A, N', N" 


Kc 


z^ 


VH 


VHS 


^VH 


9,4,1 


0.413 


1.23-1.37 


1.63-1.56 


0.50-0.56 


2.00 


16, 9, 4 


0.425 


1.16-1.22 


1.68-1.65 


0.49-0.52 


1.95 


25,16, 9 


0.430 


1.13-1.16 


1.71-1.69 


0.49-0.51 


1.92 


36,25,16 


0.433 


1.10-1.12 


1.73-1.72 


0.49-0.50 


1.91 


exact ^ 


0.441 


1 


1.875 


0.5 = 


1.875 


^Data from Ref 


s. [p7| and 


|3^] for hypercubic 


finite clusters. 






^Refs. 34,pf. 












=Ref. H. 
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TABLE II. 
The same as Table I for the d = 3 simple cubic Ising model. For hypercubic clusters one 
has A^ = L^, and for rectangular prisms A^ = L{L — 1)^ or N = L'^{L — 1) ^. 



MFRG 


N,N' 


Kc 


V 


VH 


VHS 


^VH 


2,1 


0.203 


1.50 


1.89 




2.84 


4,1 


0.207 


1.22 


2.00 




2.44 


9,4 


0.212 


1.05 


2.08 




2.18 


SBMFRG 


N, N', N" 


Kc 


V 


VH 


VHS 


vyH 


4,2,1 


0.212 


1.23-1.49 


2.00-1.91 


1.03-1.12 


2.85 


8,4,2 


0.215 


0.96-1.22 


2.15-2.00 


0.90-1.05 


2.44 


12, 8, 4 


0.201 


1.32-1.03 


1.95-2.09 


0.99-0.85 


2.15 


18,12, 8 


0.216 


1.04-1.19 


2.09-2.01 


0.95-1.03 


2.39 


27,18,12 


0.217 


0.88-0.99 


2.20-2.09 


0.86-0.97 


2.07 


other methods 


0.222 ^ 


0.63 '^ 


2.48 " 


0.8 '^ 


1.56 


'^Data from Refs. [^ 


7| and P . 










^Ref. P-4q. 












^Ref. |8^41 . 












<iRef. H. 
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TABLE III. 
The same as Table I for the critical values of the d = 2 site percolation problem. 
MFRG 



N,N' 



Pc 



■■p 



V 

Vhs 



2,1 

4,2 

9,4 

16, 9 



1/3 
0.427 
0.443 
0.473 



1.56 
1.52 
1.47 
1.43 



SBMFRG 



N, N' N" 



Pc 



p 

Vhs 



9,4,1 
16, 9, 4 



0.602 
0.585 



1.25-1.35 
1.35-1.35 



0.72-0.69 
0.66-0.65 



other methods 
^ef. [H. 
^Ref. PI. 



0.593 



1.33 
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TABLE IV. 
Critical values for the d = 2 Ising model on a square lattice according to the SBMFRG 
including the corner magnetic exponent ync ^■ 



N,N\N" 



Kr 



VH 



-VHC 



16, 9, 4 


0.469 


1.08-1.13 


1.78-1.72 


0.57-0.51 


25,16, 9 


0.460 


1.07-1.10 


1.80-1.76 


0.63-0.59 


36,25,16 


0.454 


1.06-1.07 


1.81-1.78 


0.68-0.65 



exact 



0.441 



1 



1.875 



1 



^Data from Ref. ||32| applying clusters in which the corner exponent can be defined. 

^Ref. mm. 



=Ref. |2|. 



TABLE V. Critical temperature and dynamic critical exponent for the Ising model according 
to the DMFRG. 



d=2 


d=3 


N,N' 


Kc 


z 


N,N' 


Kc 


z 


2,1 


0.347 


1.17 


2,1 


0.203 


0.97 


4,1 


0.361 


1.39 


4,2 


0.207 


1.15 


4,2 


0.370 


1.60 


8,2 


0.207 


1.32 


9,4 


0.381 


2.13 


8,4 


0.209 


1.49 




0.44P 


2.2 ±0.2^ 




0.222^ 


2.0^^ 



^Refs. 0,15 
^Ref. [^|. 

^Ref. 0-0! • 
'iRef. [^|. 
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TABLE VI. Surface dynamic critical exponent zg and critical surface coupling pc for the 
three-dimensional Ising model according to the DMFRG. 



N,N' 


zs 


Pc 


2,1 

4,2 


1.35 

1.88 


1.35 

1.48 


Monte Carlo "^ 


— 


1.52 



^Ref. ra. 



TABLE VIL Critical temperature and thermal critical exponent for the Ising model according 
to the EFRG. It is quoted in parenthesis the values from MFRG for comparison. 



d=2 


d=3 


N,N' 


Kc 


u 


N,N' 


Kc 


u 


2,1 

4,1 
4,2 


0.358 (0.347) 
0.379 (0.361) 
0.371 (0.370) 


1.39 (1.67) 
1.17 (1.45) 
1.01 (1.28) 


2,1 

4,1 
4,2 


0.206 (0.203) 
0.207 (0.207) 
0.208 (0.207) 


1.37 (1.50) 
1.24 (1.22) 
1.14 (1.27) 




0.44P 


1 ^ 




0.222"^ 


0.63^ 



^Refs. ||,|5 
^Refs. \E%m 



=Refs. [35-41 
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TABLE VIII. 
Critical temperature and thermal critical exponent for the Ising model according to the 
FSSRG. The figures for A^ = 4096 and A^' = 1024 are from Monte Carlo simulations with 
the parenthesis indicating the first uncertain digit, according to statistical fiuctuations ^. 



d=2 


d=3 


N,N' 


Kc 


V 


iV,iV' 


K, 


V 


4,2 

16,4 

4096,1024 


0.473 
0.432 

0.440(8) 


0.905 
1.053 

1.00(0) 


8,2 


0.225 


1.052 




0.4407*^ 


lb 




0.222= 


0.63 '^ 



^Data from reference [^ 
bRefs. [||,|5 



"Refs. iMi- 
'^Refs. WM^. 
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TABLE IX. 



Some 


models studied by 


means 


of the methods described in 


the text (indicated by 


"X"). 




model 




MFRG 


EFRG 


FSSRG 


NMFRG 




geometrical 




X 








Classical (static) 



pure Ising like 

random Ising like 

surface critical 

spin-glass 

random-fields 

XY and Heisenberg 

Ashkin- Teller 

ANNNI 

mixed spin 

Z(q) symmetric spin 

random D-vector 

g-state clock spin-glass 

Potts 

Blume-Capel 

lattice gas 

liquid crystal 

Tsallis Statistics 



Quantum (static) 



pure and random transverse Ising 
anisotropic Heisenberg 

random-mixture 
Dzyaloshinskii-Moryia 
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granular superconductors 



Dynamic 



Ising like 

Potts 

contact-process 

surface critical 

quantum spin 

cellular automata 
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FIGURES 
FIG. 1. (a) and (b) show schematically the smallest clusters for hypercubic lattices and (c) 

and (d) the same for the directed self avoiding random walk (DSARW) in two dimensions. Full 

circles represent sites belonging to the cluster itself and open circles the corresponding surrounding 

sites (mean field sites). 
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